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Abstract 

We analyze the local density of states (LDOS) of heterostructures consisting of al- 
ternating ferromagnet, F, and superconductor, S, layers. We consider structures of 
the SFS and SFSFSFS type, with thin nanometer scale F and S layers, within 
the ballistic regime. The spin-splitting effects of the ferromagnet and the mutual 
coupling between the S regions, yield several nontrivial stable and metastable pair 
amplitude configurations, and we find that the details of the spatial behavior of 
the pair amplitude govern the calculated electronic spectra. These are reflected in 
discernible signatures of the LDOS. The roles that the magnetic exchange energy, 
interface scattering strength, and the Fermi wavevector mismatch each have on the 
LDOS for the different allowed junction configurations, are systematically investi- 
gated. 
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1 Introduction 



The continual investigation and development of heterostructures consisting 
of layered ferromagnet (F) and superconductor (S) compounds at nanome- 
ter length scales has captured the interest of researchers in a broad range of 
disciplines. The reasons for this stem from the distinctive properties of these 
nano-sized materials, modified from those of their bulk counterparts, related 
to the effects induced by the mutual couphng that occurs when supercon- 
ductors and ferromagnets are in close electrical contact. The advancement of 
nanofabrication techniques is partly responsible for the maintained interest 
thus far, where electron-beam lithography, etching and lift-off methods, to 
name a few, have been refined to allow for the well controlled creation of rel- 
atively clean multiple layers involving F/S junctions. Further developments 
in experimental probes such as the scanning tunneling microscope (STM), 
have also progressed to the point where //eV energy resolution can be ob- 
tained at very low temperatures[l]. When operating in spectroscopy mode, 
the STM gives fine atomic-scale details of the local density of states (LDOS) 
of inhomogeneous superconducting systems. In contrast to conventional tunnel 
junctions, where transport measurements are taken, the STM gives localized 
spatial information not averaged out over the extent of the entire sample[2]. 

For composite materials comprised of multiple superconductor layers sepa- 
rated by ferromagnet materials, the presence of the magnetic F regions in- 
duces a sphtting of the Fermi surfaces of the different spin bands, and leads 
to spin-dependent Andreev[3] states. The degree of spin-polarization in the 
superconductor and the induced pairing correlations in the ferromagnet are 
all components that embody what arc known as proximity cfi'ccts. It turns 
out that due to the mutual interaction between the two contrasting materials 
and the inherent system inhomogeneity, the pair amplitude, becomes spatially 
varying[4,5], and its phase difference A0 can modulate between successive S 
layers. For certain configurations, there can be (at zero applied magnetic field) 
a phase difference of A0 = vr between S layers separated by a ferromagnet [6] . 
These are the so-called n junctions, which occur in the simplest case, in a 
3-layer SFS type structure. Junctions of this type can also occur in different 
combinations for more complicated[7.8] layered heterostructures, where the 
relative sign of the pair potential A(r) can change between adjacent S lay- 
ers. The proximity effect is then responsible for the existence of a number of 
stable junction configurations of the pair amplitude in F/S multilayers [9]: the 
oscillatory behavior of the superconducting correlations in each ferromagnet 
and the spatial profile of the pair amplitude in the superconductor are closely 
interconnected, often resulting in a nontrivial spatial dependence of the mag- 
nitude and phase of the order parameter. This consequently implies that the 
one-particle energy spectrum, or density of states (DOS), will have an equally 
interesting behavior. 
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The characteristic behavior of the DOS in structures consisting of a single 
ferromagnet and superconductor has been addressed in past work. In the dif- 
fusive, dirty-hmit case, the damped-osciUatory decay of the pair amphtude in 
the F region induces a spatial variation of the LDOS[10]. For a clean F/S 
system with rough boundaries, the tunneling DOS as a function of energy 
revealed certain inversion signatures that depended on the thickness of the 
ferromagnet film [11]. The energy gap in F/S hybrid systems becomes rapidly 
suppressed, in a manner that depends on the ratio of the F and S Fermi ener- 
gies, interfacial scattering strength, and F layer thickness[12]. The depletion 
of the gap has also been discussed in the context of a potential spin-valve 
device [13]. The spatially averaged DOS was calculated for a clean 2-D F/S 
box geometry[14], as was the LDOS for 3-D hybrid F/S nanostructures[15], 
where there was a reported filling in of the gap with increased F layer thick- 
ness. The damped oscillatory phenomena of the DOS is expected to persist 
over a broad range of material parameters [16]. For an SFS type structure, 
the proximity effect can lead to a modification of the LDOS that depends 
solely on the phase difference of the superconductors [17]. It was shown by 
using a two-band model for ferromagnetism, that the LDOS for a F/S bi- 
layer still maintained the characteristic damped oscillations, with however, an 
increased wavelength and sharper decline [18]. Certain F/S bilayers with un- 
conventional d-wave superconductors have also been investigated for peculiar 
characteristics of the DOS [19,20,21]. Thus, while the DOS and LDOS oiF/S 
bilayers have been investigated, as have also to some extent SFS junctions, 
work analyzing the LDOS of junctions consisting of more than three layers 
is still scarce. It is imperative therefore to have a theory valid for these mul- 
tilayers, since as potential device applications and experimental techniques 
become more advanced, the creation of complicated sequences of 7r-junctions 
will become increasingly realizable. 

Much of the somewhat limited experimental work on the DOS in F/S het- 
erostructures is based on transport measurements. Direct evidence of the os- 
cillatory behavior of the superconducting correlations in the ferromagnet was 
found through tunneling spectroscopy measurements that showed inversions 
in the DOS for a thin ferromagnet film [22]. Modifications to the DOS in the 
superconductor also give insight in describing the proximity effect, by pro- 
viding information on the influence of the ferromagnet on superconducting 
correlations. It was found that the DOS in the vicinity of the interface of the 
superconductor is substantially modifled from the bulk BCS result [23]. Re- 
cently, local spectroscopy revealed the LDOS for a few different ferromagnet 
thicknesses [24]: the reported inversion of the DOS was attributed to the tt 
state. Some other relevant experiments involving the ground state behavior of 
TT-junctions have shown that both the or vr state can occur, depending on 
the width dp of the F layer [25]. The origin of the damped oscillatory behavior 
reflected in DOS measurements and calculations discussed earlier is responsi- 
ble, in some cases, for the damped oscillations in the critical current Ic as a 
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function of F layer width, suggesting a to tt transition [26]. The signature in 
the characteristic Iq curves also indicates a crossover from the to tt phase in 
going from higher to lower temperatures[27]. The current phase relation was 
also measured[28], demonstrating a re-entrant Iq with temperature variation. 

Despite the existence of this experimental work, little attention has been paid 
until very recently to discussing, from a thermodynamic point of view, the 
relative stability of the different possible configurations involving and tt 
junctions in SFS . . . heterostructures. Indeed, a recent preprint [29] complains 
of the complete absence of work on thermodynamic properties of such struc- 
tures. We have alleviated this situation in very recent work[9], where we have 
found a way to correctly evaluate the condensation free energy in these nanos- 
tructures as a function of the relevant parameters. The work that we discuss 
here involves discussing the LDOS for junction configurations whose thermo- 
dynamic stability properties are[9] well understood. 

The quasiparticle amplitudes and energies that ultimately govern all observ- 
able behavior are very sensitive to the geometry of the given F/ S system. For 
F OT S layers a few nanometers thick, interference of the wavefunctions even 
over the atomic scale can produce significant contributions to the LDOS. Of- 
ten the interaction potentials and geometry vary quite rapidly over the Fermi 
wavelength, and averaging over the momentum space trajectories eliminates 
useful information on the quasiparticle dynamics. It is therefore preferred in 
the clean limit situation considered here, to implement a microscopic model 
which affords the most complete information regarding the local electronic 
structure. This requires going beyond the standard approximate equations, 
e.g., the Eilcnberger equations, which are essentially transport differential 
equations that describe superconductivity in the quasiclassical regime. With 
the progressive strides made in local spectroscopy techniques, a fully micro- 
scopic theory that allows inclusion of the details of the LDOS at the atomic 
scale is necessary. Furthermore the ground state properties of tt junctions can 
be analyzed properly only if the superconducting order parameter is calculated 
self- consistently. The theory used here satisfies these requirements. 

In this paper we calculate the LDOS for several possible multilayer structures. 
The existence of such structures as either local or global minima of the free en- 
ergy was proved in Rcf. [9]. Differing spectra are found depending on the exact 
nature of the particular junction configuration considered. By using an estab- 
lished fully microscopic method [7,30,31], we self-consistently calculate the pair 
potential, A(r), revealing several permissible types of spatial arrangements of 
the order parameter. The quasiparticle wavefunctions and associated energies 
are obtained from the microscopic Bogoliubov-de Gennes (BdG) equations. A 
numerical algorithm is then used that permits the self-consistent calculation 
of the pair potential. We find differing signatures for the LDOS that should 
be discernible in STM spectroscopy experiments. 
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2 Methods 



In this section we outline the Bogoliubov-de Gennes (BdG) equations [32] that 
are appropriate for the investigations of F/S muhilayered junctions. Since the 

BdG equations are inherently real-space based, they are especially convenient 
for the investigation of the possible spatial configurations of the pair amplitude 
and the spatially resolved LDOS. 

We consider 3-D slab-like compounds translationally invariant in the x — y 
plane, and with all spatial variations occurring in the z direction. The struc- 
tures consists of alternating superconducting, S, and ferromagnetic, F, layers, 
each of width ds and dp respectively. We call d the total thickness of the slab. 
The corresponding coupled equations for the spin-up and spin-down quasipar- 
ticle amplitudes {u}^, vD are written, 

Ef{z) + U{z) - ho{z)]ui{z) + A{z)vi{z) = enui{z) 

(la) 

E^{z) + U{z) + ho{z)]vi{z) + A{z)ui{z) = e„^;JL(z), 

(lb) 

where e± is the kinetic energy term corresponding to quasiparticles with mo- 
menta transverse to the z direction, e„ are the energy eigenvalues, A{z) is 
the pair potential, and U{z) is the scattering potential at each F/S interface, 
given by: 

(iVi-l)/2 

U{z)^ 5^ [U(i(ds + dF),z) + U((i(ds + dF)-dF),z)], (2) 

i=l 

where we take U{zi,z) = HS{z — zi), zi is the location of the interface, H 
is the scattering parameter, and A^^ represents the total number of layers 
(superconducting plus magnetic). The ferromagnetic exchange energy Hq^z) 
takes the constant value ho in the F layers, and zero elsewhere. Other relevant 
material parameters are taken into account through the variable bandwidth 
Ef{z). This is taken to be Ef{z) — Efs in the S layers, while in the F layers 
one has Ef{z) = Efm so that in these regions the up and down bandwidths 
are respectively Ep^ = Efm + ho, and Efi = Efm — ho- The dimensionless 
parameter /, defined as / = ho/ Efm, conveniently characterizes the magnets' 
strength. The ratio A = Efm/Efs = {kFM/kFsY describes the mismatch 
between Fermi wavevectors on the F and S sides, assuming parabolic bands 
with kps denoting the Fermi wave vector in the S regions. The dimensionless 
parameter Hb = mH/kFs characterizes the interfacial scattering strength. We 
consider here singlet pairing only and, consistent with this choice, we assume 
that the magnetic orientation in all magnetic layers, when more than one is 
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present, is the same. 



The nontrivial spatial dependence of the pair potential must be calculated in 
a self consistent manner by an appropriate sum over states: 

A(^) = Y: I de^ [ui{z)vi{z) + ui{z)vl{z)\ tanh(6„/2T), (3) 

\en\<WD 

where A'"(0) is the density of states (DOS) per spin of the superconductor 
in the normal state, d is the total system size in the z direction, T is the 
temperature, is the cutoff "Debye" energy of the pairing interaction, and 
g{z) is the effective coupling, which we take to be a constant g within the 
superconductor regions and zero elsewhere. The quasiparticle amplitudes, 
and involved in the sum are extracted from the solution to (1) by using 
symmetry arguments [30]. 

An appropriate choice of basis allows Eqs. (1) to be transformed into a finite 
2N X 2N dimensional matrix eigenvalue problem in wave vector space: 



H+ D 
D H- 



*n = en (4) 



where = {uli, . . . , uj^j^, v^i, . . . , w^at), are the expansion coefficients associ- 
ated with the set of orthonormal basis vectors, ul^{z) — ^j2/d Yl,q=i u\q s\n.[kqz) , 



and vlX^) = \J2/ dY^q^]^vl^qSV[i{kqz). The longitudinal momentum index kq is 
quantized via kq = q/ird, where g is a positive integer. The label n encom- 
passes a sum over the index q and the value of s±. The finite range of the 
pairing interaction cud, implies that is finite. We have calculated the matrix 
elements in Eq. (4) elsewhere [9], and for brevity, the results are suppressed 
here. 

We are mainly interested in the quantity of most experimental relevance: the 
local density of states (LDOS). The LDOS is related to the convolution of 
one-particle spectral functions for the S and F regions, and is defined as 
N{z,e) — N'^{z,e) -\- N\^{z,e), where the LDOS for each spin orientation is 
given by, 

N^iz^e) = -E|K(^)]V'(e- 6„) + K(^)]V'(e + 6„)|. (5) 
Here a [ and /'(e) = df /de is the derivative of the Fermi function. 
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3 Results 



In this section we present results for the LDOS and for the pair amphtude 
F(r) = (■0j^(r)'0|(r)), where the ijj^ are the usual annihilation operators. The 
results are given as a function of the dimensionless parameters J, Hb and A, as 
indicated. The geometrical parameters will be given in dimensionless form in 
terms of the inverse of kps'- Dg = kpsds and Dp = kpsdr- We take Dp = 10 
and ds — ^o: where is the BCS coherence length, with Ds — 100. The cutoff 
frequency is fixed to cud — O.OAEps- From these values and standard relations 
one infers g. The temperature will be fixed to T = 0.01T°, where T° is the 
bulk transition temperature of S. All results will be presented in convenient 
dimensionless form as indicated. 

Our equations (including the self consistent condition) are numerically solved 
using a previously described[7,9,30,31] iteration process. The different possible 
self consistent states are found by starting the iteration process with initial 
guesses of different types. For example, for an SFS structure one can start 
either with a "0" or a "tt" initial guess. The resulting self consistent state 
reached after iteration and convergence will be of the initial type, if such a 
self consistent state exists at least as a metastable state (a local minimum of 
the free energy) for the parameter values and geometry under consideration: 
otherwise[9] the process converges to a different configuration. 

We consider three layer SFS and seven layer SFSFSFS structures. In the 
first case, there is only one junction, and it can be in either the or the tt state. 
In the second case, there are three junctions, each of which can in principle be 
in a or TT state. Each overall junction configuration can then be characterized 
by specifying the three values, each zero or tt, corresponding to the state of 
each junction. 

3. 1 Structures 

In this subsection we present some results for the superconducting correlations 
in SFS and SFSFSFS structures, as described by the pair amplitude, F{z). 
We do this to illustrate the possible self-consistent pair amplitude configura- 
tions that arise as parameters such as the interface scattering strength Hb-, the 
Fermi wavevector mismatch A, and the magnetic polarization strength / vary. 
This short discussion will establish what the different structures encountered 
look like, some of their stability properties, and our notation. The presentation 
is most clearly and efficiently illustrated by means of 3-D plots. 

To begin, we display, in Fig. 1, the pair amplitude (normalized to Fq, its bulk 
superconductor value) for a three layer SFS system, as a function of position 
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Fig. 1. The pair amplitude F{Z), normalized to its bulk value Fq, for a three 
layer SFS structure, plotted as a function of Z = kpsz and of the dimensionless 
magnetic polarization /, at A = 1 and Hb = 0. The Z = point is at the center 
of the structure. Panel (a) corresponds to self consistent results obtained (see text) 
with an initial guess where the junction is of the "0" type and panel (b) with a "vr" 
type. In the latter case, the solution found is of the vr type except at very small 
/ and in the neighborhood of / ~ 0.4, but in the former case (panel (a)), a more 
complicated behavior occurs, as discussed in the text. 

and of the magnetic exchange parameter / in their entire ranges, without a 
barrier or mismatch (A = 1 and Hb = 0). The top panel shows the results of 
attempting to find a solution of the type by starting the iteration process 
with an assumed piecewise constant of that form. Whether this attempt suc- 
ceeds or not is reflected, of course, in the self consistent results plotted. Careful 
examination of the two panels reveals a rather intricate situation: one can see 
in the bottom panel that a solution of the vr type exists nearly in the entire 
/ range, the exception being at very small /, where the effect of magnetism 
becomes very weak and, as one would expect, only the state solution is 
found. This requires small values of J, J < 0.1 however. One can see that near 
this small value, as the state transitions to a vr configuration, the amplitude 
of F{z) is very small throughout the sample, which should be reflected by a 
sharp dip in the transition temperature. A brief instability of the vr state also 
occurs in a small interval around / ~ 0.4, where the self consistent solution 



Fig. 2. The normalized pair amplitude F{Z) for a SFSFSFS seven layer structure, 
plotted as a function of Z and of the barrier height parameter Hb- Here / = 0.2 
and A = 1. Panels (a) (b) (c) and (d) correspond, respectively, to attempts to find 
solutions of the 000, vrvrvr, ttOtt and OttO types (see text). The resulting structure, 
however, it is not necessarily of the sought type, since some of the structures are 
unstable in part of the Hb range, as explained in the text. For clarity, the direction 
of increasing Hb is not the same in every panel. 

is of the type. This instability is reflected in Fig. 7 of Ref . [9] . On the other 
hand (panel (a)), one can see that there is a region of / values, 0.12 < / < 0.3, 
where a state does not exist. The attempt to find it fails, and the solution 
converges, after a very large number of iterations, to the same self consistent 
TT state that is found, in the same / range, in panel (b). 

We now illustrate the behavior of F{z) in seven layer SFSFSFS structures, 
which can be viewed as consisting of three adjacent SFS junctions. In the 
notation described above we denote as "000" the junction configuration when 
adjacent S layers always have the same sign of A(z), and as "tttttt" when 
successive S layers alternate in sign. There are up to a trivial reversal, two 
more symmetric states: one in which A{z) has the same sign in the first two 
S layers, while in the last two it has the opposite sign, (this is labeled as the 
"OvrO" configuration), and the other corresponding to the two outer S layers 
having the same sign for A{z), opposite to that of the two inner S layers: 
these are referred to as "vrOvr" junction configurations in our notation. We 
will mainly focus our study on these symmetric configurations. For the first 
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Fig. 3. The normalized pair amplitude F{Z) for a SFSFSFS seven layer structure, 
as in the previous Figure. The top and bottom panels correspond to initial attempts 
to find the asymmetric vrOO and vrvrO states respectively. Clearly such attempts only 
succeed when the barrier is high and the proximity effects weak (see text). 

example we consider the effect of the barrier thickness, as determined by the 
parameter Hb- This is done in Fig. 2. In this figure we have taken / = 0.2 
and A = 1 (no mismatch). The geometrical parameters and temperature are 
the same as in Fig. 1. The panels (a), (b), (c) and (d) correspond (in this 
order) to attempts to find solutions of the 000, tttttt, vrOvr, and OvrO type. 
One can see that in all cases each attempt succeeds when Hb is sufficiently 
large, (the direction of increasing Hb has been reversed in two of the panels 
so that this feature is not concealed by the perspective). This makes sense, 
as at larger barrier heights the layers become more isolated from each other 
and the proximity effects must eventually become insignificant. The relative 
phase of A(z) at each S layer should then become progressively irrelevant. 
For the same reason, there is a clear increasing trend with increasing Hb for 
the absolute values of F{Z) in the middle of the S layers. At smaller values 
of Hb, a more careful examination of the panels reveals a rather intricate 
situation: two of the junction configurations, tttttt and ttOtt (panels (b) and 
(c) respectively), exist in the entire range of Hb, while attempts to find a 000 
structure (panel (a)) result in a ttOtt configuration for small barrier heights 
Hb ^ 0.48, and similarly, attempts to find a OttO structure result in a self 
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Fig. 4. The normalized pair amplitude F{Z) for a SFSFSFS seven layer structure, 
plotted as in Fig. 2 and with the same panel arrangements, but as a function of the 
mismatch parameter A at zero barrier. Note that the direction of increasing A is, 
for clarity, not the same in all the panels. 

consistent tttttt configuration for Hb ^ 0.44. This is in perfect agreement with 
the stability analysis of Ref . [9]. Both of these "transitions" consist simply of 
the flipping of the two outer junctions, with the inner one remaining in the 
same configuration. 

For completeness, we show in Fig. 3, results for the two possible asymmetric 
states that can exist (vrOO and vrvrO). The parameter values are identical to the 
symmetric case exhibited in the previous figure. The top panel (a) corresponds 
to the self-consistent profile when an assumed form for the pair potential cor- 
responds to a vrOO junction, while the bottom panel is for a nnO initial guess. 
It is evident that, for reasons specified above, at large enough values of the 
scattering parameter Hb both asymmetric states are stable self-consistent so- 
lutions. We find that the vrOO state is stable for Hb ^ 0.52, while the ttttO 
junction is more robust, retaining its asymmetric configuration over a larger 
range of Hb values, Hb ^ 0.35. However, for weaker barriers, where the prox- 
imity effects are important, these asymmetric states are not stable, becoming 
TTTTTT and vrOvr respectively. 

In Fig. 4 we consider the influence of the Fermi wavevector mismatch, charac- 
terized by the parameter A = Efm/Efs- This figure is completely analogous 



11 



Fig. 5. The normalized DOS (see text) for a SFS trilayer, plotted as a function of 
the energy (in units of the bulk zero temperature gap), and of the exchange energy 
parameter / = ho/Eps- The panel arrangement is the same as in Fig. 1: therefore 
the plots correspond to and vr self consistent states as in that figure. 

to Fig. 2, except for the substitution of A for Hb-, which is set to zero (no 
barrier). We see that for sufficiently large mismatch (small A), all four junc- 
tion configurations exist. This is because increasing the mismatch does, in 
effect, isolate the superconducting regions from each other just as efficiently 
as increasing the barrier does. This can be roughly understood qualitatively 
from simple quantum mechanical considerations. For smaller mismatch, only 
the two solutions of the tttttt and ttOtt type exist: again 000 turns into ttOtt 
(for A > 0.48) while OvrO turns into tttttt at A > 0.57, these results being in 
agreement [9] with those of detailed stability studies. 



3.2 DOS 



We now present some of our results for the LDOS, which is the main focus of 
our paper. Although the quantity F{z) gives useful information regarding the 
superconducting correlations in the structure, it is the LDOS which is experi- 
mentally measured. The results given are in all cases for the quantity N„{z, e) 
defined in Eq. (5), summed over a, and integrated over a distance of one co- 
herence length (or, equivalently, the thickness of one superconducting layer) 
from the edge of the sample. We normalize our results to the corresponding 
value for the normal state DOS of the superconducting material, and the en- 
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Fig. 6. The normalized DOS for a SFS trilayer, plotted as a function of the energy, 
and of the mismatch parameter A, at / = 0.2 and Hb = 0. The panel arrangement 
again corresponds to that in Fig. 1. In the bottom panel, the results plotted are for 
a self consistent vr structure, while in the top panel they are for a structure for 
A < 0.64 and a vr structure at smaller mismatch, when the structure is not stable. 

ergies to the bulk value of the gap at zero temperature, Aq. The results shown 
are for the three and seven layer geometries discussed at the beginning of this 
section. Material parameters not otherwise mentioned are set to the "default" 
values of Hb = 0, A = 1 and / = 0.2 

Figure 5 characterizes the sensitivity of the LDOS, in a SFS junction, to varia- 
tions in the dimensionless exchange field parameter, I. The panel arrangement 
is as in Fig. 1. The top panel, therefore, shows the resultant DOS when the 
structure is predominately in the state configuration, the exception being a 
small range of I as discussed in conjunction with Fig. 1, where the state is 
unstable. The bottom panel shows the results for the vr junction case, except 
at very small I. The limiting case of an SNS junction (/ = 0) displays (as 
can be seen in either panel) subgap Andreev bound states that quickly decay 
as / increases. Upon varying J, the number of subgap states oscillates as a 
function of /, in a manner that is dependent on whether the stable state is 

or TT. 

In Fig. 6 we show some other results for the same three layer system. The figure 
shows the normalized local density of states, as specified above, as a function 
of energy and mismatch A, at fixed / = 0.2. The panel arrangement is again as 
in Fig. 1: the bottom panel represents the results for a n configuration, while 
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Fig. 7. The normalized DOS for a SFS trilayer, plotted and displayed as in Fig. 6, 
and under the same conditions, but as a function of barrier height rather than 
mismatch. The results shown in the bottom panel correspond to a vr structure, 
while those the top panel they are for a structure for Hb > 0.29 and a vr structure 
at smaller barrier height values, when the structure is not stable (see text). 

the top panel shows the results for the configuration at larger mismatch, 
A < 0.64, while for smaller mismatch, where the configuration does not [9] 
exist, it repeats the results of the vr configuration. The change in the results 
refiecting the onset of the configuration instability at A = 0.64 can be 
discerned in panel (a). At larger mismatch (smaller A, the back portion of the 
figures), where both states can coexist, one can clearly see in the figure the 
difference between the two configurations. This is helped because in general 
when there is more mismatch the proximity effect weakens, the gap opens and 
larger peaks form which progressively become more BCS-like as the mismatch 
increases (A decreases) to A = 0.2. The way this occurs, too, is not the same 
for the and the n configurations. 

In the next figure. Fig. 7, we repeat the previous study but in terms of the 
barrier height, rather than the mismatch. Again, we plot in each panel the 
self consistent results obtained with initial guesses of the and vr types, which 
means that all regions of the (b) panel correspond to self consistent vr states 
whereas those in the top panel correspond to configurations, except for 
Hb < 0.29 where no self consistent state exists[9]. We then plot instead 
results for the same vr state as in the bottom panel, which is reached by 
iteration of either an initial or vr guess. Again, the point where the transition 
occurs can be seen as a clear discontinuity in panel (a). The difference between 
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Fig. 8. The normalized DOS for a SFS trilayer, plotted as a function of the energy. 
Results for both and vr self consistent states are given, as indicated. These are 
slices from Figs. 6 and 7. In the top panel, the DOS is shown at Hb = for two 
different values of the mismatch parameter, A = 0.2, and 0.67, the latter being a 
case for which the state is nearly unstable (see text). The bottom panel shows the 
DOS profile in the absence of mismatch (A = 1), but with the interface scattering 
parameter Hb taking on the two values shown, chosen on similar criteria as the A 
values (see text). 

the DOS for the two states where they coexist is in this case even more obvious 
than in the previous one, with the BCS like peaks being considerably higher 
for the 7T configuration. 

It is also illustrative to display more clearly some of the features of the results 
by isolating some selected slices of the 3-D plots as 2-D line drawings. This 
we do in Fig. 8. In the top panel of this figure we present results for an 
SFS trilayer, for two contrasting values of the mismatch parameter A. Results 
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Fig. 9. The normalized DOS for the SFSFSFS structure of Figs. 2 and 4, with 
the same panel arrangements as in those figures, plotted as a function of the energy 
and of the mismatch parameter A. 

labeled as "0" and "vr" are for the case where the self consistent states plotted 
are of these types. The and n state curves corresponding to A = 0.67, 
where (as shown above and in Ref. [9]) the state is barely metastable, have 
clearly distinct signatures, with a smaller gap opening for the state, and 
consequently more subgap quasiparticle states. Consistent with this, when 
there is little mismatch one finds that the pair amplitude is relatively large 
in the F layer. In agreement with what was seen in Fig. 6, this progression 
takes a different form for the and vr states, and this is reflected in this 
panel, as can see by careful comparison of the curves. In the bottom panel 
we demonstrate the effect of the barrier height. The results are displayed as 
in the top panel, but with the dimensionless height Hb taking the place of 
the mismatch parameter. This figure should be viewed in conjunction with 
Fig. 7. One of the values of Hb chosen {Hb = 0.32) is again such that the 
state is [9] barely metastable, while for the other value the and vr states have 
similar condensation energies. For Hb close to 0.32 there is a marked contrast 
between the two plots, with the gap clearly opening wider and containing more 
structure for the vr state. At larger Hb the gap becomes larger in both cases, 
with the BCS-like peaks increasing in height. Thus, as the barrier becomes 
larger, one is dealing with nearly independent superconducting slabs and the 
plots become more similar. The largest difference therefore occurs, as for the 
mismatch, at the intermediate values more likely to be found experimentally. 

Now we turn to the DOS results for the seven layer structure. In Fig. 9 we 
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Fig. 10. The normalized DOS for the SFSFSFS structure of Figs. 2, 4 and 9, 
plotted as in Fig. 9 except as a function of barrier height Hb instead of A. 

present results for the normalized local DOS, defined as in the three layer case, 
as a function of energy and mismatch parameter. The panel arrangement in 
this figure corresponds to that in Fig. 4, with however a lower bound of A = 0.3 
in the mismatch parameter, in order to better display the more intricate LDOS 
structure. The results shown are for the self consistent state reached by starting 
from each of the four symmetric junction configurations. For larger mismatch, 
which corresponds to the back portion of the plots, all results are therefore 
similar, as the LDOS reflects the independent nature of the isolated slabs (the 
phase is irrelevant in this limit). For small mismatch (the ranges are indicated 
in the discussion of Fig. 4) the results of panel (a) coincide with those of panel 
(c) and those of panel (d) with those of (b), as can be seen by comparing 
the front portions of these panel pairs. The transitions can be noted in the 
plots. Examination of the results shows clear differences between the four 
configurations. As in the three layer case, the peaks are more prominent at 
larger mismatch, and the "gap" is wider. At small mismatch, where only two 
configurations are stable, the gaps fill up and the peaks become very small for 
both cases, namely tttttt (panel (b) or large A region of panel (d)), and vrOvr 
(panel (c) or front region of panel (a)). 

A similar situation occurs as a function of the barrier parameter Hb as shown 
in Fig. 10, which corresponds to the same conditions as Fig. 2, with the same 
conventions and panel arrangements. In this case, however, the BCS peaks are 
(at larger barrier heights) more prominent for the tttttt and ttOtt configurations 
since the S layers are not yet completely decoupled at the largest Hb shown. 
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which (again for clarity) is smaller than that used in Fig. 2. The peaks weaken 
and the gap narrows and fills up at small barrier heights, as the proximity 
efi:ect strengthens. The prominence or weakening of the peaks is associated[9] 
with an increase or decrease of the condensation energy of the corresponding 
configuration, as one would intuitively expect. For small barrier heights we 
note also that the two states that exist in this case (tttttt and ttOtt) have again 
clearly different signatures, as can be seen by comparing the front (small Hb) 
portions of diagrams (a) and (c) with those of (b) and (d). A similar situation 
occurs at small mismatch (larger A), in Fig. 9. 



4 Conclusions 

We have seen in this paper that the different junction configurations that 
compete in stability for different ranges of the relevant parameters, yield rather 
characteristic signatures in the local DOS. Thus, our results illustrate that the 
complex and nontrivial behavior of the pair amplitude in clean layered F/S 
nanostructures is refiected in the LDOS. The numerical method used here 
is based upon self consistent solution of the BdG equations. One advantage 
of this method is that it properly takes into account the atomic-scale and 
quantum interference effects that arc likely to be important in these systems. 
For the SFS and SFSFSFS structures investigated here, one can have up 
to two and six independent locally stable configurations, respectively; as the 
number of layers increases, so too does the number of possible configurations. 
Quantitatively, we have shown that the larger the scattering strength at the 
interfaces, or the mismatch between Fermi wavevectors of the F and S regions, 
the greater the number of stable states, since the proximity effects weaken. The 
various pair amplitude profiles however, exhibited distinct spatial forms in each 
case. We confirmed[7,9,33] that it is essential to calculate the pair potential 
self consistently, since [9] as the results shown here demonstrate, multilayer 
configurations can transition from one state to another abruptly, as one varies 
any one of the physical parameters: the magnetic exchange energy /, interface 
scattering strength Hb, and the Fermi wavevector mismatch A. 

The diverse behavior observed in the pair amplitude is found to be reflected 
in the calculated energy spectra. When several self-consistent pair amplitude 
configurations coexist, each representing a local or global minimum of the free 
energy, they yield specific signatures in the local DOS, which are susceptible to 
measurement in local tunneling spectroscopy experiments. For a trilayer SFS 
junction, variations in the exchange field parameter / lead to oscillations in 
the number of subgap states, and there are similar changes in the seven layer 
spectrum. We also investigated the LDOS profiles by separately varying Hb 
and A. In the limiting cases of large Hb or small A the LDOS takes its typical 
BCS-like form, with few subgap quasiparticle states. The LDOS in general 
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however, has characteristic peaks and dips that depend on the symmetry of 
the ground state, and on Hb, I or A. 

It would be of obvious interest to examine whether or not the configuration 
changes that occur as the material parameters are varied can happen also as 
a function of temperature. This question is now under investigation. 
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